Stochastic Population Dynamics Driven by 

Levy Noise 



Jianhai Bao and Chenggui Yuan* 
Department of Mathematics, 
Swansea University, Swansea SA2 8PP, UK 



Abstract 

This paper considers stochastic population dynamics driven by Levy noise. The 
contributions of this paper lie in that (a) Using Khasminskii-Mao theorem, we show 
that the stochastic differential equation associated with the model has a unique global 
positive solution; (b) Applying an exponential martingale inequality with jumps, we 
discuss the asymptotic pathwise estimation of such model. 
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1 Introduction 

Stochastic population dynamics perturbed by Brownian motion has been studied exten- 
sively by many authors. There are a great amount of literature on this topic. In [6j Mao, 
Marion and Renshaw investigate stochastic n- dimensional Lotka-Volterra system 



dX(t) = diag(X 1 (t),--- ,X n (t))[(b + AX(t))dt + aX(t)dW(t)}, (1.1) 

where 

X = (Xl, • • • , X n ) , b — (bl, • ■ ■ ,b n ) , A = (ajj) n xn, Cr = (Cij)nxn, 

and reveal the important role that the environmental noise can suppress a potential pop- 
ulation explosion; Mao, Sabanis and Renshaw [7J further discuss the asymptotic behaviour 
of population process determined by Eq. (II. ip . Then the techniques developed in [61 [7J 
have been applied successfully to study stochastic delay population dynamics in [21 [9] ; 
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functional Kolmogorov-type systems [HI [12]; hybrid competitive Lotka-Volterra models in 

[131 mug. 

The population may suffer sudden environmental shocks, e.g., earthquakes, hurricanes, 
epidemics, etc. However, stochastic extension of population process described by Eq. (11. ip 
cannot explain the phenomena above. To explain these phenomena, introducing a jump 
process into underlying population dynamics is one of the important methods. To our 
knowledge there is few systematic work so far in which the noise source is a jump process. 
The work presented here is to take some steps in this direction, building extensively on the 
existing results mentioned above for the Brownian motion case, is also a sequel to that of 
[3], where competitive Lotka-Volterra population dynamics with jumps 

dY(t) =Y(t)\(a(t) -b(t)Y(tj)dt + a(t)dW(t) + [ j(t,u)N(dt,du) 



' JY 

is investigated, and the explicit solution, sample Lyapunov exponent and invariant measure 
are also addressed. 

We focus in this paper on stochastic population dynamics ( II .ip that is further perturbed 
by Levy noise, that is, 

dX(t) = diagpTi(t), • • • , X n (t)) \{b + AX(t))dt + aX(t)dW(t) 

(1.2) 

+ / H(X(t-),u)N(dt,du) 



Here b,A,a are defined as in model (11.11) . W(t) is a scalar Brownian motion defined on 
the probability space {f^J 7 , P} with the filtration {J r t}t>o satisfying the usual condition, 
X(t~) := lim s ^ t X (s) , N(dt, du) is a real-valued Poisson counting measure with characteristic 
measure A on measurable subset Y of [0, oo) with A(Y) < oo, N(dt, du) := N(dt, du) — 
X(du)dt, and H : M. n x Y — > W 1 . Throughout this paper, we further assume that W is 
independent of N. 

In reference to the existing results in the literature, our contributions are as follows: 

• We use jump diffusion to model the evolutions of population dynamics when they suffer 
sudden environmental shocks; 

• Using Khasminskii-Mao theorem, we show that the Stochastic Differential Equation 
(SDE) associated with the model has a unique global positive solution; 

• Applying an exponential martingale inequality with jumps, together with the standard 
Borel-Cantelli lemma, we discuss the asymptotic pathwise estimation of such model. 



2 Global Positive Solutions 



Since X denotes the population sizes of the n interacting species, it is natural to require 
the solution of Eq. (II. 2p not only to be positive but also not to explode in a finite time. 
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Therefore, in this section we intend to show that Eq. (11.21) has a unique global positive 
solution under some conditions. Since the coefficients don't satisfy linear growth condition 
or weak coercivity condition, even they satisfy local Lipschitz condition, the solutions of Eq. 
(11.21) may explode in a finite time. Khasminskii [4J Theorem 4.1, p85] and Mao [5] gave the 
Lyapunov function argument, which is a powerful test for nonexplosion of solutions without 
linear growth contion and is referred as Khasminskii-Mao theorem. In what follows, we shall 
also apply Khasminskii-Mao approaches to show that Eq. (11 .2p has a unique global positive 
solution X(t),t > 0. In this section we will show: 

• Jump processes can suppresses the explosion; 

• Brownian motion can also suppresses the explosion to our new model, which is similar 
to that of [6]. 

In what follows let K > be a generic constant whose values may vary for its different 
appearances. To show the main result let us recall the following facts. 

Consider 1-dimensional SDE with jumps 

dX(t) =F(X(t))dt + G(X(t))dW(t)+ [ <$>(X{r),u)N(dt,du), t>0 (2.1) 

Jy 

with initial condition X(0) = xq G R, where W(t) is a real-valued Brownian motion, F, G : 
l^R,$:lx¥^l. 

The following conclusion is given by (see [101 Lemma 2]. 

Lemma 2.1. Let F(0) = G(0) = $(0,w) = for u G Y and F, G, $ satisfy local Lipschitz 
condition. Set 

. i r /, \x + $(x, u)\\ 2 w , . 

J(x) := J ^ln^ J - n J \(du),Xj£0. 

Assume that 

sup J(x) < oo for each m > 1. (2.2) 

0<|x|<m 

Then for x ^ 

F(X(t, x)^0 and X(t~, x) ^ for any t > 0) = 1, 
where X(t, x) denotes the solution of Eq. (12.11) starting from x at time t = 0. 

For convenience of reference, we recall some fundamental inequalities stated as a lemma. 
Lemma 2.2. 

x r < l + r(x- 1), x>0, l>r>0, (2.3) 

n 

n (i-f)A0| x | P < J- x P< n (i-|)vO| x | P? Vp>0) i6 k; (2.4) 
i=i 

where K" := {x G 1™ : X; > 0, 1 < i < n}, and 

lna; < x — 1, x > 0. (2.5) 
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2.1 Explosive Suppression by Jump Processes 

For jump- diffusion coefficient we assume that 
(HI) For any m > 1, x G R n , «6Y and i — 1, • • • , n 

H i {x,u)>-1, fl i (0,«) = 0, (2.6) 
sup / (In |1 + Hi(x,u)\) 2 \(du) < oo, (2.7) 

;|z|<m Jy 



0< 

and for each > there exists constant L k > such that 



/ \H(x,u) - H(y,u)\ 2 X(du) < L k \x - y\ 
Jy 



{21 



whenever x,y G R™ with |x| V \y\ < k. 



Theorem 2.1. Let assumption (HI) hold. Assume further that for p e (0,1) there exist 
constants 5 > 0, a > 2 such that for x G M™, z = 1, . . . n, 

J^p) := / [(1 + fl- i (x,u)) p - 1 -pfli(x,u)] \{du) < -5\x\ a + o(|x| a ), (2.9) 
Jy 

where o(|x| a )/|x| a — > as \x\ — > oo. Then, for any initial condition x G M™, Eq. fll.2p has a 
unique global solution X(t) G for any t > almost surely. 

Before we prove the theorem, we give an example such that condition f)2.9p holds. 

Example 2.1. For i — 1, ■ • ■ , n, let 7, > and J Y (1 V 7 f (w))A(d«) < 00. Assume that 

Hi(x,u) := ii(u)H(\x\), xeR n ,ueY, 

where if(|a;|) is a polynomial of degree a > 2 with positive leading coefficient. Then by a 
straightforward computation we have for some S > 



JKx,p)= / [(l + ji(u)H(\x\)y-l-pji(u)H(\x\)]\(du) 
Jy 

< -S\x\ a + o(\x\ a ). 
Therefore condition (12.91) holds. 

Proof of Theorem 12.11 By (12.81) . for arbitrary initial value x G R" there is a unique local 
solution X(t) for t G [0, r e ), where r e is the explosion time. By Eq. (11.21) the ith component 
Xi(t) of X(t) admits the form 

n n 

dXiit) = Xiit) (bi + OijXj(t)\dt + Y VijXj{t)dW{t) 

3=1 ' 3=1 (2.10) 

+ / H i (X{r),u)N(dt,du)] 



Note that for any t G [0, r e ) 

i 1 n 

(&i + ( s ) - 2 ( 

i=i i=i 



+ / (ln(l + iZitX'(s), it)) - u))X(du) )ds 

+ / V<7yX i (s)dW(s)+ / / ln(l + i^(X(0,w))iV(ds,du)). 
Jo - =1 Jo Jy ^ 

This, together with x G R", yields that > for any t G [0, r e ). On the other hand, 

due to 11277)) . for Eq. (127IU1) condition (E2D holds. Then Lemma O gives X 4 (t) > for any 
t G [0, T e ) since x G . Next we show r e = oo a.s. Let fco > be sufficiently large such that 
\x\ < k . For each k > k define a stopping time 

r fc :=inf{tG [0,T e ) : \X(t)\ > k}. 

Clearly, is increasing as k f oo. Set r ra := lim^oo r^, whence < r e a.s., then it is 
sufficient to check = oo a.s. Introduce a Lyapunov function for any p G (0, 1) 



i=i 

Let T > be arbitrary. For any < t < A T, the Ito formula yields 

n n 

dV{X(t)) = CV(X{t))dt + p^X p {t) ^aijXj(t)dW(t) 

i=l j=l 

n „ 

+ V / [(l + H i (X(r),u)) p -l]N(dt,du)Xf(t), 

i=l JY 



(2.11) 



(2.12) 



where for x G 



i=i 



6^ ~t~ ^ ^ (ZijXj 



1 — p 



3=1 



v3'=l 



CijXj 



p 

1; 



By inequality f )2M| 
and, thanks to (12. 9p 



+ V / [(l + ^(x, M ))P-l-j9^(x, M )]A(rf 

i=l JY 

K"i(x,p) + K 2 (x,p). 

Ki(x,p) < K\x\ 2+P + o(\x\ 2+p ), 
K 2 (x,p) < -5\x\ a+p + o(\x\ a+p ). 



(2.13) 



U )X; 
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Thus, for a > 2 

CV(x) < K for any x £ R". (2.14) 

Define for each u > 

//(«) := influx), |x| > u}. 
Thanks to inequality f )2.4p . it is easy to see that 

lim fi(u) = oo. (2-15) 

it— >oo 

Then we obtain from (I2.14p that for some constant K > 

fi(k)¥(r k < T) < E(V(X(r k ))I {Tk < T} ) < EV(X(r k A T)) < K. 
Recalling (I2.15P and letting k — )• oo yields 

P(roo < T) = 0. 

Since T is arbitrary, we must have 

P(r oc = oo) = 1 

and Eq. (11.21) admits a unique global solution X(t) e R™ on t > 0. 

Remark 2.1. In [6], under the condition 

(H2) an > if 1 < % < n while a {j > if i ^ j. 

Mao, Marion and Renshaw reveal the important fact that Brownian motion noise can sup- 
press a potential population explosion. Theorem 12.11 shows that Levy noise can also play the 
same role, without any conditions being imposed on the diffusion coefficient a. 

As for population dynamics, in general, the following Lyapunov function for p e (0, 1) 

n 

U(x) := J2i xP i ~ 1 -Plnxi], xeRl, (2.16) 
i=i 

is constructed to show that the SDE associated to the model admits a unique global positive 
solution, see, e.g., [SI [131 O US]- In what follows, under suitable conditions we can also 
show that Eq. (jl.2p has a unique global positive solution through the Lyapunov function 
defined by (|fl6j) . 

Theorem 2.2. Suppose that assumptions (12.61) . (12. 8ft and ( 12. 9ft hold. Assume further that 
there exist constants (5 6 (0, a] and z/ > such that 

/ [Hi{x,u) -]n(l + Hi(x,u))]X(du) < u\xf + o(\xf) (2.17) 
Jy 

for i — 1, • • • , n and x £ 1R™ . Then, for any initial condition a; £ IR™ , Eq. (11.21) has a unique 
global solution X(t) £ for any t > almost surely. 
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Proof. Since the argument is similar to that of [6, Theorem 2.1], we here only sketch the 
proof to point out the variation from the Brownian motion case. Let k$ G N be sufficiently 
large such that every component of x is contained in the interval (^, ko). For each k > k 
define a stopping time 

T fc := inf jt G [0, r e ) : Xi(t) (^r, k^j for some i — 1, • • • , n|, 

where r e is the explosion time. In the sequel, we show := lim^oo = oo a.s. Let T > 
be arbitrary. For any < t < r k A T, applying Ito's formula, we obtain 



CU(x):=pJ2 

i=i 



+ V / [(1 + Hi(x, u)) p - 1 - P Hi(x, u)]X(du)x p 
i=i - /y 

+ pV / [fii(ar,«) -ln(l + i?i(x,u))] A(d«) 



:= h{x,p) + h{x,p) + h{x,p). 
By flUD, note that 



< fsT|a;| 2+p + o(| 



Also, due to (EHD and f l2T7D 

J 2 (x,p) + h(x,p) < -5\x\ a+p + npu\xf + o(\x\ a+p ) 
Then for p G (0, 1) and G (0, a] 

CU(x) <K, x G M™. 

Define for each u > 1 



:= inf < £7(x) : ij > m or < - for some z = 1, • • • , n 
I u 

Due to the property of function h(x) := x — 1 — In x, x > 0, we see that 

lim /i(a;) = oo and lim h(x) = oo 

ccfoo xi.0 

and hence 

lim fi{u) = oo. 

u— ¥oo 

The proof is then complete by repeating the procedure of [6, Theorem 2.1]. 



(2.18) 



(2.19) 
□ 
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2.2 Explosive Suppression by Brownian Motion 



In this subsection, we further show that Brownian motion can also suppresses the explo- 
sion to our model under condition (H2), but weaker conditions imposed on jump- diffusion 
coefficient. 

Theorem 2.3. Assume that assumptions (HI) and (H2) hold. Then, for any initial con- 
dition x G R™, Eq. fll.2p has a unique global solution X(t) G R™ for any t > almost 
surely. 

Proof. Since the proof is very similar to that of Theorem I2.1[ we here only give an outline of 
the argument. In (I2.13p . note from (H2) and inequality (12 .4p that for p G (0, 1) and x G R™ 

Ki(x,p) < _ P^LlPh 2 min au \x\ 2+p + o(|a;| 2+p ). (2.20) 

2 l<i<n 

On the other hand, by inequality (I2.3p . for any p G (0, 1) and x G R™ 

/ [(1 + Hi(x, u)f - 1 - pHfa u)} \{du) < 0, 
Jy 

hence we have 

K 2 (x,p) < 0, p G (0, 1) and x G R". 

Thus 

CV{x) < K for any x G K". 

The conclusion then follows by carrying out the procedure of Theorem 12.11 □ 

Applying Lyapunov function U(x) in (I2.16p . under suitable conditions we can still guar- 
antee that Eq. (II. 2p admits a unique global positive solution. 

Theorem 2.4. Let conditions (12. 6p . (12. 8p . (H2) hold and assume further that condition 
(12.171) holds with (3 G (0, 2]. Then, for any initial condition x G K." , Eq. (I1.2p has a unique 
global solution X(t) G for any t > almost surely. 

Proof. The argument is similar to that of Theorem 12.21 In Eq. (I2.18P , by inequality (12. 3p 

I 2 {x,p) < 0, p G (0, 1) and x G K". 
What's more, condition (I2TT7]) leads to, for f3 G (0, 2] and x G R+, 

h{ x iP) < npv\x\^ + o(|x|^). 
Combining (I2.20p . for p G (0, 1) and (3 G (0, 2], we can conclude that for some K > 

CU(x) <K, x G R™. 

Then the proof can be done by carrying out the procedure of Theorem 12.21 □ 
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Remark 2.2. By constructing two different Lyapunov functions V(x) and U(x) defined by 
(12.111) and (12.161) . respectively, under different conditions, namely, (12. 7\i and (I2.17p . we show 
that Eq. (11.21) has a unique global positive solution. Comparing conditions (12.71) and (12.171) . 
we see that condition (12. 7p is a local one while (I2.17P has a growth restriction on jump 
diffusion H in whole space. Therefore, the argument developed in Theorem 12.31 is easier to 
verify than that of Theorem 12.41 

3 Asymptotic Moment Properties 

In the last section, under suitable conditions we have shown that Eq. (I1.2p admits a 
unique global positive solution. From the biological point of view, the nonexplosion property 
and positive solution in a population dynamical system are often not good enough while the 
moment properties are more desired. In this section we shall show that the pth moment 
with p G (0, 1) and the average in time of the moment of the solution to Eq. (II. 2p are both 
bounded. 

Theorem 3.1. Under conditions of Theorem 12.11 (or Theorem 12.21 or Theorem 12.3} or 
Theorem 12.41) . for any p G (0, 1) and some K > 

limsupE|X(t)| p < K. (3.1) 

t— >oo 

Proof. We only give a sketch of proof under conditions of Theorem 12.11 due to the similarities 
of arguments. For any p G (0, 1) let V be defined by (12. lip . For any < k define a stopping 
time 

a k := inf{t > : \X(t)\ > k}. 

By the argument of Theorem 12.11 we have <Tfc < oo and 07. t 00 a - s - as k — > oo. Applying 
the Ito formula yields 

E(e tAu "V(X(t A <7 fc ))) = V(£(0)) + E / e s [V(X(s)) + CV(X(s))}ds, (3.2) 

Jo 

where CV is defined by (12.131) . Since the leading terms of CV are — <5|x| cH " ? ' with a > 2 we 
can deduce that there exists a constant K > such that V{x) + £V(x) < K. Hence 

E(e'l/(X(t))) < ff(l + e*) 

and the desired conclusion follows from inequality (12 ,4p . □ 

Going through the arguments of Theorem 12.11 (or Theorem 12. 2\ or Theorem 12.31 or 
Theorem 12.4)1 . we can also derive that the average in time of the moment of the solution to 
Eq. (II .2p is bounded. 

Theorem 3.2. (i). Under conditions of Theorem 12. II or Theorem 12 .2\ for any p G (0, 1) and 
any initial value x G there exists constant K > such that 

1 r l 

limsup- / E\X(s)\ p+2 ds < K. 
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(ii). Under conditions of Theorem 12.31 or Theorem 12. 4} for any p e (0, 1) and any initial 



value x e W! , there exists constant K > such that 



1 r* 

limsup- / E|X(s)| p+2 ds < AT. 

t-s>oo I Jo 



Proof. Since the proofs of (i) and (ii) are very similar, we here only prove (ii) under the 
conditions of Theorem 12.31 The argument is motivated by that of [7J Theorem 2}. By (I2.12p 
and (I2.13P we have 



/t / n n n i ~\ n \ 

(pE^+pEE x ^ oo - E "S^oo ) 
v j = i j = i 7= i j = i / 



V(X(t)) < F(X(0)) + 

+ Mi(t) + M 2 (t), 

where Mi(t),M 2 (t) are two local martingales. Noting that the polynomial 

n n n / 1 \ n 

v E 6 * + p E E ~ ^4 — E a ii x ^ +p 

i=l i=l j'=l i=l 

has a upper bound .K" (dependent on p), therefore 

v(x(t)) + £ < v ( x (°)) + Kt + + M ^(t). 

i=l 

Taking the expectation, and then dividing by £ on both sides, we obtain the result. □ 
By taking another different Lyapunov function, we will have the following theorem. 

Theorem 3.3. Under the conditions of Theorem 12. 1[ let p T = (pi, . . .p n ) be positive num- 
bers such that pi + . . . + p n < 1, and assume there exist constants j5\ and fa G (0, a) such 
that 



nr =1 (i + mix, «))» - x)a + «)) Pl 

-' L i= i 



A(d«) < fa\xf 2 + o(\xf 2 ), (3.3) 



where constant a > 2 was given in (12.91) . Then 

E(IT7 =1 Xf (t)) < oo, for all t > 0. 

Remark 3.1. For n = 1, condition (13.31) must be true. Moreover, Example 12.11 also demon- 
strates that condition (13. 3p holds in some cases. 

Proof of Theorem MM Define a C 2 -function V : M™ -»■ M + by 

V(x) :=nr =1 xf. 
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Compute 



1. 



CV(x) = V(x)p T (b + Ax) — -V (x)x T a 1 [diag(pi, . . . ,p n ) — pp T ]ax 



nr =1 (i + Hi(x,u)) pi - 1 - y^PiHi 



x, u) 



X(du) 



V(x)p T (b + Ax) — -V(x)x T a T [di&g(pi, . . . ,p n ) — pp T \ox 



[l + Hi&u^-l-piHi&u) 



X(du) 



+ V(x) J2 [ 

f r n 
+ V(x) / n™ =1 (l + Hi(x, u))' H - V(l + Hi(x, u)) Pi + n - 1 

L i=1 



A(du). (3.4) 



Noting conditions (12. 9 j) and f 13 . 3 1) . we derive that there exist positive constants Ci and C*2 
such that 

< V r (x)(C 1 -C 2 |x| a ). 
For each fc > \x\ define a stopping time 

T k := inf{t > : \X{t)\ > k}. 

By Theorem 12.11 r k < oo and — )■ oo as fc — >■ oo almost surely. Using the Ito formula we 
obtain 

"tAr fc 



EV(X(t A r k )) = V(X ) + E / CV(X(s))dt 

Jo 



< V{X ) + CiE / V(X(s))ds. 
Jo 

Hence applying the well-known Gronwall inequality and letting k — > oo gives 

EV(X(t)) < V(X )e c ^, 
and the required assertion follows. 



4 Asymptotic Pathwise Estimation 

In the last section we have discussed how the solutions vary in in probability or in 
moment. In this section we examine pathwise properties of the solutions. To discuss the 
pathwise properties of Eq. (11.21) . we cite the following exponential martingale inequality 
with jumps, , e.g., [H Theorem 5.2.9, p291]. 

Lemma 4.1. Assume that g : [0, oo) — y R and h : [0, oo) x¥->R are both predictable 
J-" r adapted processes such that for any T > 



T pT 

2 



\g(t)\ dt < oo a.s. and 



o Jy 



\h(t, u)\ 2 X(du)dt < oo a.s. 
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Then for any constants a, (3 > 



F\ sup 

0<t<T 



a 



g (s)dW{s)-^ I \g{s)\ 2 ds + 



h(s, u)N(ds, du) 



JY 



[ e «h(.,«) _ j _ a/l ( S) M )]A(dw)cfe > /3 1 < e" Q/3 . 
a Jo Jy j j 

Theorem 4.1. Let conditions of Theorem 12.11 hold. Assume further that there exists con- 
stant 9 G (0, a] such that 



/ [(In Q (x, u)f + Q(x,u))X(du) < K\x\ e + o(| 
Jy 

where, for p G (0, 1) and x G M+ 

n n 

Q(s, it) := £(1 + «))^ / £ a* 



(4.1) 



(4.2) 



i=i 



There exists > 0, independent of initial value x G K", such that the solution X(t),t > 0, 
of Eq. (11. 2p has the property 



ln(|X(t)|) r ^ 
limsup u ^ JU < K, a.s. 
t->oo In t 



(4.3) 



Remark 4.1. By inequality ( 12. 4ft . it is easy to see that condition ( 14. 1]) holds for ifj in 
Example 12.11 

In what follows we complete the argument of Theorem 14.11 

Proof of Theorem 14.11 Note by Theorem 12.11 that Eq. ( II .2p has a unique global positive 
solution for any initial value x G E". Let V(x), Ki(x,p) be defined by ( 12. lip and ( 12 . 13j) . 
respectively, and Z(x) := YH=i P x ^i Y^=i < 7 ij x j f° r x e ^+ anc ^ P e (0' -0- Applying the 
Ito formula to e* In V(x) yields 

e'lnViXit)) = lnV(x) + j\ s [\nV(X(s)) + * - ^ 2 (X(s)) 



+ / lnQ(X( S ), M ) 



P 



53xf(s)£T i (X a ,u))A(dM) 



(4.4) 



</¥ 



+ / e s Z(X(s))dW{s) + / / e s lnQ(X(s-),M)iV(ds,du) 



By virtue of Lemma I4.1[ for any a, (3,T > we have 



F\ u) : sup 

0<t<T 
+ 



e s Z(X(s))cW(s) 



a: 



e 2s Z 2 (X(s))ds 



e s kiQ{X{s-),u)N{ds,du) 
~aJ J (Q ae3 ( X ( s )^ u )- 1 ~ aeSln Q( X ( s )^ u )) X ( du ) ds >P}< e 



o Jy 

t 



-a/3 
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Choose T = k,a = ee k and (3 = 2e lnk , where k G N, < e < \, in the above equation 



2' 



Since Y2T=i k 2 < oo, we can deduce from the Borel-Cantalli lemma that there exists an 
Q C Q with P(f2o) = 1 such that for any uj G Qq we can find an integer ko(u) > such that 



f e s Z(X(s))dW(s) + [ [ e s \nQ(X(s-) lU )N(ds,du) 
Jo Jo Jy 



< + ^ f e 2s Z 2 {X{s))ds (4.5) 



2 Jo 

+ e^Io I { QeeS ^ X ^ s ^ u ^~ 1 ~ eeS ~ klnQ ^ X ^ s ^ u ^) X ^ ds 
whenever < t < k and k > k (u). Hence, for any uj G ffo, < t < k and k > ko(ou) 



\nV(X(t)) <e- t \nV(x) + 

1 - 1 , , l~e w2 



+ / e s - 1 



\nV(X(s)) + vix(s)) K i( X (s),P) - —Z\X{s)) 



ds 



eT*^ (\nQ(X(s),u) - yfj^fi J2xf(s)H l (X(s),u))\(du)ds 

+ _ Tfc I I (Q eeS ~ k {X{s),u) -l-ee s ~ k In Q{X{s),u))\{du)ds 
ee 'o Jy ^ ' 



:= Ji(t) + J 2 (t) + J 3 (t) + J 4 (t). 
For any x G R™ and -u G Y, compute that 

n 

In Q(x, u) - -7— V a^iZ^x, u) 

n 

= logQfx, «) - Q(x, u) + 1 + w) - — - — xZHJx, u) - 1 

Km ^—^ 

(4-6) 

< Q(a;, u) - — — a%Hi(x, u) - 1 

^ ' i=l 

1 " 

where in the second step we used the inequality (12.51) . In the light of a Taylor's series 
expansion, for e G (0, |], x G R™ , w G Y and s < fc 

Q eeS_fc (x,M) = 1 + ee s ~ k \nQ(x,u) H (lnQ(x, u)) 2 Q^(x, u), 

where £ lies between and e. Thus 

J 4 (t)= I I ^—^ (\nQ(x,u)) 2 Q^x,u)X(du)ds. 



JY 
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Note that for any u G fi , t < k and k > ko(uj) 

n2s-t-k 
(In Q(x, u)) 2 Qt(x, u)\{du)ds 
_<Q(x,«)<l 2 

n ee 2s~t~k 
(In Q(x, u)) 2 Q^(x, u)X(du)ds 
J(x,u)>l 2 

= : r 1 (t) + r 2 (t). 

For < Q(x, u) < 1 and < £ < e < |, we have Q^(x, u) < 1. Hence 

n2« f 
^— (In Q(x, w)) 2 A(dw)ds 
„<QO,u)<l 2 
« 2s t k 

< — (\nQ(x,u)) 2 X(du)ds. 

Jo Jy 2 

On the other hand, recalling the fundamental inequality 

In a; < 4(x^ — 1) for x > 1, 
and observing Q^(x,u) < Q^(x,u) for Q(x,u) > 1 and < £ < e < |, we have 

£ & 

^— Q(x,u)X(du)ds 

4(x,u)>i 2 

Q(x, u)X(du)ds. 

2 

Consequently, for any u <E Q , t < k and k > k (cu) 

Mt)<- I [ e s "'[(lnQ(x,M)) 2 + lQQ(x,u)}X(du)ds. 
2 Jo Jy 

Then, for any u £ Qq, < t < k and k > ko(u>), we have 

Mt) + Mt) + J 4 (t) < jf e-* [in V(X(s)) + —L—^^^.p) - ^Z 2 (X(s)) 

+ ^ / [(lng(X( S ), M )) 2 + 16Q(X( S ), M )]A(rfM)lc/ S 

e s -'[M 1 (X( S ),p) + M 2 (X(s),p) + M 3 (X(s),p)}ds. 

By inequality (12. 4p and inequality (12. 5p . for x G R+ and p G (0, 1) we obtain 

Mt(x,p) < K\x\ 2 + o(\x\ 2 ). 

14 



o 



Moreover, again by inequality (12. 4p . together with (I2.9P and (14. ip . we can deduce that for 
x e M+ and a > 2 

M 2 (x,p) + M 3 {x,p) < -K\x\ a + o(\x\ a ). 

Noting that Mi(x,p) + M 2 (x,p) + M 3 (x,p) is bounded by a polynomial with the negative 
leading coefficient, we arrive at 

J 2 (t) + J 3 (t) + J 4 (t) < K [ e'-'ds = K{1 - e~*). 

Jo 

Thus, for any u G fi , < t < k and k > ko(uj) 

In V(X(t)) < e~* In V(x) + — ^ + K. 

In particular, for u £ Qq, k — 1 <t < k and k > fco(w) + 1, we have 

M\X(tm lnV(X(t)) 1_ + + 2e>n* 

In -f In I . . 1 1 L J 



\nt ~ hit ~ln(A;-l) L w J eln(A;-l)' 

This implies that 

ln(|X(t)|) 2e 
hmsup < — . 

t^oo m t pt 

The desired assertion then follows by letting e t f • 

Noting the limit limt^oo ^ = 0, we can easily deduce from Theorem 14. II that the sample 
Lyapunov exponent of Eq. (11. 2p is less or equal to zero, which is stated as the following 
corollary. 



Corollary 4.1. Under conditions of Theorem 14. 1[ 

limsup Vl v JU < 0, a.s. 

t— >oo t 

Remark 4.2. Our theories developed can also be applied to discuss stochastic functional 
Kolmogorov-type population dynamics with jumps and switching- diffusion ecosystems with 
jumps, respectively, which will be reported in separated papers. 
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